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We study the problem of difTusing particles which coalesce upon contact. With the aid of a 
non-perturbative renormalization group, we first analyze the dynamics emerging below the critical 
dimension two, where strong fluctuations imply anomalously slow decay. Above two dimensions, the 
long-time, low-density behavior is known to conform with the law of mass action. For this case, we 
establish an exact mapping between the physics at the microscopic scale (lattice structure, particle 
shape and size) and the macroscopic decay rate in the law of mass action. In addition, we identify 
a term violating this classical law. It originates in long-range and many-particle fluctuations and is 
a simple, universal function of the macroscopic decay rate. 
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I. INTRODUCTION 

The standard for the mathematical treatment of chem- 
ical reaction kinetics is provided by the law of mass action 
(LMA), stating that the rate of an elementary reaction 
is proportional to the product of the densities of its par- 
ticipants [HE]- However, unless the environment is con- 
tinually well-stirred, the interplay between reaction and 
dispersion implies fluctuations in space and time, which 
tend to promote retarded kinetics [SJH]. Indeed, typically 
there is a critical dimension below which ineffective diffu- 
sive mixing leads to a qualitative change in the speed of 
the reaction, and one speaks of anomalous kinetics [T]|S]- 
[7]. In contrast, above the critical dimension, one usually 
expects fluctuation only to alter the decay kinetics quan- 
titively, by influencing the reaction rates f2j. Yet, rigor- 
ous analysis often only supports the validity of the LMA 
for asymptotically long times and low densities, whereas 
analytical treatment away from this regime is rare [HI [S] ■ 

Employing a non-perturbative renormalization group 
(NPRG) approach, in this work we explore a bimolecu- 
lar reaction scheme where particles clot upon contact, 
A + A A, and explore space by diffusion. In the 
spirit of investigating idealized and simplified models to 
gain insight into microscopically much more complex pro- 
cesses, this scheme was introduced by Smoluchowski to 
help understand the physics behind the coagulation of 
gold particles suspended in an electrolyte [TUHT^ . This 
type of reactions is of interest for a variety of fields of 
natural science, with applications to the dynamics of 
aerosols [T3], the binding of proteins [5], the decay of 
defects in solids [7], or the kinetics of chemical reac- 
tions PE]. 

From perturbative renormalization group analysis it is 
known that, below the critical dimension two, the reac- 
tion kinetics displays critical behavior with anomalously 
slow decay of the density |15l [T5] . We recover this re- 
sult, and demonstrate that for the one-dimensional case 
the NPRG allows to obtain a substantially enhanced 
quantitative result as compared to the perturbative ap- 
proach. The NPRG is also suited for the analysis above 



the critical dimension. We can thus show that for three- 
dimensional systems the LMA is violated by a relatively 
strong, non-analytic term in the density. The effects be- 
come more pronounced as the particle density increases 
and are attributed to long-range fluctuations and many- 
particle correlations. The corresponding additional term 
is a universal function of the non-universal, macroscopic 
decay rate (the proportionality factor in the LMA). To 
calculate the latter, we construct a mapping between the 
physics at the microscopic scale, defined by the shape and 
size of the particles and the structure of the lattice, and 
an effective macroscopic description. To the best of our 
knowledge this is one of the very few instances where the 
renormalization group idea of mapping physics from the 
micro- to the macro-scale can be accomplished exactly. 

The rest of this article is organized as follows. In Sec- 
tion [n] we give a brief introduction to the NPRG for- 
malism. Starting from the microscopic definition of the 
process, it allows to calculate the effective action F, by 
systematically integrating over the degrees of freedom, 
going from high to low momenta and frequencies. From 
F the macroscopic behavior, such as the kinetic equa- 
tion governing the time evolution of the density, can be 
deduced. After these preliminaries, we proceed in Sec- 
tion |TTT] to discuss the coagulation process below the crit- 
ical dimension. We demonstrate that the specific prop- 
erties of the problem at hand imply significant simplifi- 
cations of the flow equation. This enables us to derive an 
adequate numerical description of the anomalous decay 
kinetics, characterized by universal parameters. Finally, 
in Section |IV| we turn to the treatment of the process 
above the critical dimension in a broader and more de- 
tailed discussion than is given in our previous paper on 
this topic [17j . 

II. NPRG AND REACTION-DIFFUSION 
SYSTEMS 

The NPRG based on the effective action F is a suitable 
framework for the study of fluctuation effects in reaction- 
diffusion systems. Originally developed for the treatment 
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of equilibrium physics |18H21| . it has been adapted re- 
cently to non-equilibrium systems |22ff5i] . In the con- 
text of reaction-diffusion, it has led to new insights in 
the critical properties of branching and annihilating ran- 
dom walks and in the generalized voter class [25j. 

For completeness and to set the notation we provide 
in this section a short introduction to the NPRG for- 
malism. We specialize the treatment to the coagulation 
process and put particular emphasis on the derivation 
of the kinetic equation. For authoritative reviews and 
more detailed accounts on the general topic the reader is 
referred to [MMl M 



A. The Field-Theoretic Action 

In a first step, the process is mapped onto a field the- 
ory by a Fock space formalism, an approach devised by 
several authors |S1 E51 - I5T] , for a review see • In our 
case of diffusing particles which coagulate upon contact, 
the action of the path integral can be split up as 



(1) 



Here the fields (f> — 0(x,i) e IR and (/> = (l){x,t) e ilR 
are related to particle annihilation and creation opera- 
tors, respectively, and are variables of the position x on 
the lattice and of time t. There is a term for the time 
evolution 

where the factor Z = 1 {it may change along the renor- 
malization group flow, although for the coagulation pro- 
cess this not the CclSGj ctS discussed below). Furthermore, 
we have a diffusion term Se[(f>, (j)] (a particle can hop from 
site X to site y with rate h, when x and y are adjacent 
sites, indicated by brackets < x, y > in the following 
sum), which becomes 



and we set the lattice spacing equal to one, a = 1, to 
define the length scale. In this case the hopping rate h 
equals the diffusion constant D. To define the scale of 
time, we set, unless otherwise stated, D — 1. 

Below the critical dimension, the structure of the lat- 
tice does not infiuence the long-time, low-density behav- 
ior [ini . For the purposes of Section III it is therefore ad- 



equate to perform the continuum limit, where e(q) 
i.e. Eq. ^ is replaced with 

J q,LJ 

dxdt^(x,t)VV(x,t) . 



Finally in the action (jlj there is a reaction term 
S\[4>, (/)] for coagulation A + A ^ A, which reads 

J2 [dt A(y - x) [0(y, t) + l] 0(x, t)0(y, t)0(x, t) . 
x,y 

In order to be able to study interactions with a finite 
range, we have introduced here the reaction kernel A(z), 
which defines the shape and size of the particles. Thus, 
if there is a particle at site x, a particle at site y is anni- 
hilated with rate A(y — x). 

Again, below the critical dimension the microscopic 
properties are irrelevant for the universal behavior and it 
suffices to consider local interactions 



Sx 



^(x,t)0(x,i)2 



with A := '^i'^)- 



B. The Wetterich Equation 

In analogy to equilibrium physics, for reaction- 
diffusion systems one can define a generating functional 



Zh J2 /di [0(x,i)-^(y,O] [(/'(x,t)-(^(y,t)] Z[J,J]=J V^V(bexp ~S[^,(j)\+^j 



<x,y> 



dt {J(j) + Jcf) 



in position space (again Z — \). In Fourier space, due to 
translational invariance, it is diagonal. 



z 



e(q)</'(-q, -w)0(q,a;) 



(2) 



where f :— f^^, and f :— ( -tA runs over the first 
Brillouin zone of volume Vb- Unless otherwise stated, 
in this article, we consider a hypercubic lattice, where 
Vb — (^) and the dispersion relation reads [33] 



e(q)=4/.^sin2(^) , 



from which observables such as the temporal evolution 
of the particle density can be derived. The fields J, J are 
introduced in order that the associated functional 

W[J,J] = lnZ[J, J] 

genererates the connected Green's functions by func- 
tional derivation in the fields J, J at J J = 0. We 
remark that the field J induces particle input — > ^, as 
long as J is positive. 

The degrees of freedom are integrated systematically, 
going from fast frequencies and short wavelengths at the 
microscopic scale A (related to the typical reciprocal 
length scale of the interaction) to slow frequencies and 
long wavelengths. This is implemented by introducing a 
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"mass term" ASk to the action suppressing fluctuations 
of these modes, 



Zk[J,J] = / I?0X'</)exp 



SI 



J2 ldt{J(f> + J^) 



where the mass term reads 

ASk[^,(j>]^ i?fc(q)(^(-q, -w)0(q,w) , 

J q,aj 

with a cutoff function Rk that is independent of the fre- 
quency w, a convenient choice for reaction-diffusion pro- 
cesses [12] . Specifically, we take [33J [31] 

i?fc(q)-(fc'-e(q)) e{k^ - e{q)) . 



(3) 

Instead of calculating the renormalization group flow 
of the functional Wk[J, J] = lnZfe[J, J] [351 |3B], for the 
approach employed in this article one considers the so 
called effective average action Ffe. Denoting the expecta- 
tion values 

:= (0) = SWk/5J , V (0) - SWk/SJ , (4) 
at J = J = 0, it is deflned by 

rfc[V^, ^] = Y1 (^^ + '^^) - ^'^fcl'/'' - ^fcl"^' ^] • 

(5) 

Notice that this is just the Legendre transform of Wk, 
up to the term ASk, which is added for mathematical 
convenience and vanishes as fc ^ 0. 

At the microscopic scale A the cutoff function ASk is 
supposed to be large, so that it freezes the deviations 
from the expectation values ip, and and renders the 
functional trivial. The initial condition thus becomes 
(see ^24j for a careful discussion) 

Only modes with 9 ^ A are integrated out along the 
renormalization group flow. Therefore, if we employ the 
continuum limit in the action 5 (as in Section III I, a fl- 



nite A sets the ultraviolet cutoff. If otherwise we keep the 
discrete lattice structure, q only runs over the first Bril- 
louin zone, intrinsically enforcing an ultraviolet cutoff. 
Similarly, we have an implicit ultraviolet cutoff for finite- 
size objects, such as balls with radius R, where modes 
with q ^ R~^ are suppressed. In these cases there is no 
need for an explicit cutoff and we can set A = cx). 

The central equation of the NPRG formalism is the 
Wetterich equation. This flow equation connects the ef- 
fective average action at the microscopic scale T^^a — S, 
where fluctuations are neglected, with the effective action 
r = Tk=o where all degrees of freedom are integrated. It 
can be expressed as 



dkTk[^,^]^^Tr 



Ri 



(6) 



" (2) 

where and Rk denote the 2x2 matrices of the second 
functional derivatives of Tk and of the mass term AS'^, 
respectively, and Tr denotes the trace. 



C. The Derivative Expansion 

In general, the Wetterich equation cannot be solved 
exactly. The derivative expansion is an approximation 
based on the fact that infrared singularities are sup- 
pressed by the cutoff function Rk and therefore the ef- 
fective average action rfc['0,'0] is analytic so long as the 
scale fc > [55]. It is performed by expanding the func- 
tional Tklip, ip] in orders of the temporal derivative dt and 
the spatial derivative V, truncating Tk[ipi 4'] at a certain 
order. 

For the following expansions, let us take the continuum 
limit of the action as the initial condition. A common 
truncation of the effective average action is the "leading 
order" approximation 

rfe[^,V^]= Jd''xdt[Uki^,^p) + Zk^p{dt-DkV^)^] . 

(7) 

This approximation is popular because it already allows 
to determine the anomalous dimension 77 = — fc9fclnZfc 
and the dynamic exponent z ^ 2 + kdk In Dk ■ How- 
ever, we show below that for the coagulation process this 
ansatz is equivalent to a lower order in the approxima- 
tion, the so called "local potential" approximation, 

Tki^,^]^ [ d^x dt [Uk (V^, V-) + ^ {dt - V^) 1^] , (8) 



where only the renormalization of the local potential Uk 
is taken into account. The reason is that, as proven in 
Section |III| for our particular process Zk and Dk are 
not affected by the renormalization group flow, i.e. Zk = 
Dk — I for all k. Thus, for the coagulation process rj and 
z agree with the mean-fleld exponents and 2, respec- 
tively. 

Let us consider homogeneous flelds 'ip{x,t) = 
'0, ip{x,t) = Ip. The relation between the average ef- 
fective action and the local potential then reads 

Tk[^,^]^VT-Uk{^,tP), 

with the asymptotically large volumes of space and time 
V ^ Jd-^x^ {2tt)'^6{p = 0) and T = Jdt ^ 2nS{uj = 0), 
respectively. With the cutoff function (Isl (where e(q) = 
, since we consider the continuum limit) , from the Wet- 
terich equation ^ one can deduce the flow equation for 
the effective average potential. 



dkUk 



(9) 



Here Vd = 



Vd denotes the volume of the d- 



dimensional unit sphere, and 



(m,n) 



the mth 
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and nth derivative of Uk with respect to ij) and respec- 
tively. In our case the initial condition reads 

D. Dimensionless Flow Equation 

In order to be able to study critical behavior, we need 
to resolve the fixed points of the flow. As verifled in the 
next section, this is realized if we introduce the dimen- 
sionless coordinates 

x = fc^"'^x, < = fc~^i , q=fcq, uj = k^Cj , (10) 

and the renornialized dimensionless fields 

7^(x,t) =x(x,?), i^M = k\{i,t). (11) 

Also, we introduce the renormalization time t = ln(A;/A). 
The proper dimensionless and renormalized form of the 
local potential obeys 

c/fe(v;,v) = fc''+V(x,x)- 

The cutoff function ([s]) now takes the form Rk{q) = 
q^r{(f') with T-independent 

rm - - l) © (1 - ■ (12) 

The flow for the dimensionless potential follows from 
Eq. (|9|, 



drUr = -{d+2)Ur + dxu'^r'^^ 

V, {ui'^'^ + 1 



41^1) + 1 -u?'"^u^^^'^ 



, (13) 



with the initial condition 

Ur=o{x,X) = Ar=OX^X^ + Ar=OXX^ , 

where Xr=o = A'^^^A is the rescaled coagulation rate. 
E. The Kinetic Equation 



Here we have allowed for particle input with rate 
J{^, t) ^ 0. The first equation is fulfilled by -0 = 0. Iden- 
tifying ip with the particle density p, the second equation 
then determines the kinetics of the process. 

The extremal principle is the macroscopic analog of the 
classical field equations, dS/Sip ~ {), 5S / 5{p — J, which 
give the mean-field equation for the particle density. No- 
tice that in contrast to the "microscopic" action S, in 
the effective action F, fluctuations in the particle num- 
bers and correlations in space and time are taken into 
account. 

When fc > 0, the effective average action can be 
written in an expansion in the fields ip, if), and multiple 
derivatives in space and time thereof [2^]. For the pur- 
poses of this article, it is sufficient to restrict the time 
evolution to homogeneous fields, i.e. -(/'(x, i) = il^it). This 
discards derivatives in space. Neglecting also derivatives 
in time except for the term ipdiipi which is already present 
in the initial action, the general form of the effective av- 
erage action reads 

with the local potential Uk, which can be defined for 
constant fields -0, -0 by F^ [0, "0] = VT Uk {ip,ip), and with 
the volumes of space and time V = (for a hypercubic 
lattice with unit lattice spacing), T = Jdt. The extremal 
principle (14) then yields 







<5F[0,p] 



5V(x,t) 



d^Uk=o{il^, p)\^=Q + dtp . 



With the non- equilibrium force 

F(p) = d^Uk=o{i^,p)\^^o 
this gives the kinetic equation 

dtp{t) = -F{p{t)) . 



(15) 



The non-equilibrium force F can be determined in 
the simulations by introducing homogeneous particle in- 
put with rate J. The extremal principle then implies 
dtp = —F{p) + J, such that for stationary states the 
non-equilibrium force equals the input rate, F{p) = J. 



As compared to the alternative non-perturbative ap- 
proach of [551 [55] an advantage in calculating the effec- 
tive average action F[-0, 0] instead of W[J, J] — lnZ[J, J] 
is that it allows direct access to the observable 0j, the av- 
erage density of the particles. At fc = 0, when all degrees 
of freedom are integrated out, the cutoff term vanishes, 
AiSfc^o = 0, and the effective average action equals the ef- 
fective action, F^^o — As a consequence of Eqs. (4pl, 
the equation of motion is obtained by the " extremal prin- 
ciple" 



5F/^0 = O, 5V/5^ = J at0 = O, V-^p. (14) 



III. RENORMALIZATION BELOW THE 
CRITICAL DIMENSION 

Although the coagulation process does not display a 
phase transition, its long-time approach to a vacant sys- 
tem can be described within the framework of critical 
phenomena. Indeed, the process is suitable for treat- 
ment with the perturbative renormalization group ap- 
proach, as was demonstrated in the pioneering works of 
Peliti [15 and Lee (TB]. Peliti established that the process 
displays an upper critical dimension and that its value is 
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dc = 2, confirming predictions based on heuristic argu- 
ments and computer simulations O El [37] . Below this 
critical dimension the density decay p ~ Ad{Dt)^^ (for 
some dimension dependent amplitude Ad and diffusion 
constant D) is significantly retarded by fluctuations as 
compared to the classical behavior p ~ A{Dt)~^. Lee 
showed that the decay amplitude Ad is amenable to per- 
turbative renormalization group analysis near the critical 
dimension, 



Ad 




21n(87r)-5 
16tt 



0{e) 



if d = rfc = 2 , 
if d < dc, 



where e — 2—d is the difference between the upper critical 
dimension dc — 2 and the dimension d. 

In this section, we first study the mathematical prop- 
erties of the effective average action Tk , exploiting special 
properties of the coagulation process. Similar to restric- 
tions due to certain symmetries in, say, magnetic models, 
we show that in the Taylor expansion of the effective av- 
erage potential Uk many terms are not generated along 
the renormalization group flow. This can most conve- 
niently be seen upon representing the flow by one-loop 
Feynman diagrams. We then proceed to exploit the rel- 
ative simplicity of the flow equations for the study of 
the coagulation process in one dimension. The calcula- 
tion cannot be carried out exactly and we need to recur 
to an approximation scheme, as introduced in the pre- 
vious section, in order to reduce the complexity of the 
flow equation. However, the coagulation process permits 
us to go to a relatively high order of the approximation, 
therefore promising accurate results. Indeed, they com- 
pare well to the exact solution for one dimension. We 
also extend our analysis to general dimension d < 2 and 
thus reproduce results from perturbative calculations for 
small e. 



A. The One-Loop Expansion and Restrictions on 
the Flow Equation 

Symmetries which are obeyed by the effective aver- 
age action usually play a central role in the dis- 
cussion of the critical behavior of a system. Exam- 
ples include the 0{N) symmetries in equilibrium statisti- 
cal mechanics |26j . and, in non-equilibrium systems, the 
KPZ-symmetries |38} [39| . the time-reversal symmetries 
of "model A" |24l HO] , or the so called "rapidity symme- 
try" of the contact process [55', 1531. Similarly, one can 
make statements on the mathematical properties of the 
effective average action Tk for the coagulation process. 
They do not come as an invariance with respect to a 
symmetry transformation of the fields ip and ip (or x and 
x)- Rather, they become apparent upon expanding the 
effective average action Tk , which is rendered analytic by 
the infrared cutoff as long as the scale fc > |55] . 

Our goal in this section is to calculate the di- 
mensionless effective average potential Ut(XiX) — 



k "^Ukiipjip)- For a general reaction-diffusion process, 
it can be expressed as a power series |24j 



E 



1 



7,lnl 



For the concrete calculations later in this section, we em- 
ploy the flow equation ( 13 1. To observe the specific prop- 



erties of the potential Ut, however, it is favorable to revert 
to the full Wetterich equation It can be recast in a 
form that allows for a diagrammatic analysis and thus 
reveals the mathematical structure more immediately. 



dk -Tr 
2 



In 



Rk 



(16) 



where the derivative dk ■— dkRk ■ dn^, acts only on the 
fc-dependence of the cutoff function Rk ■ 



The functional Vk on the right hand side in Eq. ( 16 ) is 



known from perturbative analysis as the generator of one- 
loop Feynman diagrams |26| |4T|. Therefore, the renor- 
malization group flow of the (m, n)-point vertex functions 

F^™'"' (obtained by taking m and n functional deriva- 
tives of Tk with respect to ijj and ip, respectively, at zero 
fields = ip = 0) can be represented by the one-loop 
Feynman diagrams for the (m, n)-vertex. In this way, 

we also obtain the fiow of the coefficients gi™'"-*, which 



are just proportional to F 



(7n,n) 
k (pi,LJi;...;p„ 



(where the 



functional derivatives are taken with respect to ip(pi,(jji), 



*e {1, 



i} and 'ip{pj,ijj), j e {m + 1, 



,m + n 



}) at 



zero momenta and frequencies. 

Let us consider the fiow of a general (m, n)-point ver- 
tex function, determined by the sum over the correspond- 
ing one-loop diagrams. Due to causality, the propagator 
only connects earlier to later vertices. This fact drasti- 
cally restricts the number of possible Feynman diagrams. 
Since in the initial action the number n of incoming legs is 
larger than the number m of outgoing legs for all non-zero 
vertex functions, it is impossible to construct one- loop di- 
agrams with m > n from these vertices. Therefore the 
fiow of the vertex functions with m > n is zero. Similarly, 
the minimum number of incoming legs in the Feynman 
diagrams, which is n = 2, and the minimum number of 
outgoing legs, which is to = 1, is inherited to all scales k. 

Now, from the fact that the number of legs cannot in- 
crease along the time arrow, it is readily deduced that 
the fiow of T^k'^'' vanishes, since there is no diagram of 
the form of Fig. [T]^a). As a consequence, similar to the 
absence of propagator renormalization in perturbative 
renormalization |l5j , the factors Dk and Zk are constant 



Dk 



Zk 



Therefore, the leading order approximation Q and the 
local potential approximation ([8| are equivalent for the 
coagulation process. 

Not only can we rule out certain vertex functions, 
but we can also make statements on the functional 
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tiine 

< 




Figure 1. Important one- loop diagrams determining the flow 
of the vertex functions r^™'"^ , where m and n are the num- 
ber of outgoing and incoming legs, respectively. The fact the 
number of legs can only decrease along the time arrow puts 
signiflcant restrictions on the possible diagrams. Diagram (a) 
is not created in the flow equation because the (2,l)-vertex 
function is always zero. It follows that the propagator is not 
renormalized. Diagram (b) implies a term linear in F^^'^' to 
the flow of F^^'"*'. Similarly there arise linear contributions to 
the flow of general F^™'"' (except for m — n = 2, where there 
is a quadratic term in F^^'^'). Since diagram (c) is composed 
only of (2,2)-vertices and propagators, the flow for F^^''^' can 
be solved without knowledge of further vertex functions. 




' 

1 2 3 4 5 6 
incoming legs n 

Figure 2. Illustration of the order for calculating the non-zero 
flxed-point coefficients In each step the result is in- 

dependent of the coefficients to follow. This can be shown 
by looking at all possible one-loop diagrams which determine 
the flow of the vertex function F^™'"' with n incoming and m 
outgoing legs. For the physically relevant line with m = 1 all 
the coefficients y*''" ' with n' < n must be known. In con- 
trast, the diagonal elements gf*("'") only depend on -* 
with n' < n. 



dependence of r^,'"'"-' on other vertex functions. We 
first note that the one-loop diagrams for the flow of 
an (to, n)-vertex function evidently must not contain 
(to', n')-vertices with n' — m! > n — m (the number of 
legs at one vertex in the diagram cannot decrease by more 
than the overall decrease n — m in the number of legs). 
Moreover, for m < n (except for m ~ n — 2), the flow 



dkT 



(m,n) 



is linear in F 



(m,n) _ 



For m < n the corresponding 



vertex has only m outgoing lines, which cannot connect 
again to a vertex with n incoming lines. For m = n > 2 
one loop would not suffice to include a second (m, n)- 
vertex in the diagram. In general, one-loop diagrams 
for the flow of the vertex function r^™'"' which contain 
one (to, n)-vertex can only involve exactly one additional 
vertex, which must be a (2, 2)-vertex. (For the purpose 
of illustration, in Fig. [ijb), a one-loop diagram to the 
(l,3)-vertex is shown.) Therefore, this linear term does 
not vanish. 

Let us apply our flndings to the dimensionless potential 
Mt-(XjX) with the Taylor expansion 



Ur(X,X) 



E : 

m>l,n>2.m<n 



1 



Within the local potential approximation we have that 



-,(m,n) 



T(m,n) 



(E^P^ 



fe (pi,LJi;...;p„+„,w„ + „) ~ fe (O,0,...,0,0) 9t 

0, ^jWi — 0). Thus, we can calculate the coefficients 
~*{m,n) q£ flxed-point potential u* within the local 
potential approximation step by step as follows (also see 



Fig. bl): We start with ^*(2^2)^ which is easily obtained 

f2 2) f2 2) I 1/ 

because drCjT ' depends only on Qt ' , cf. Fig. Ill c). We 
then turn to g*^^''^\ whose flow drgi^''^'' is a function of 
(jr ' and (jr ■ Assuming that we know the flxed-point 
values of ^*(™'") for all to < n we can go on to treat 
~*{m,n) successively for to = n, n — 1, . . . , 1. In each step 
one simply needs to solve the linear equation 

= ci (to, n) -f C2 (to, n) • , 

given some ci(to, n) and C2(to, n) ^ 0. More pre- 
cisely, since C2{m,n) — c' {m,n)g'^''^^ gi"^'^"^ with positive 

(2 2) 

c'{m,n) and, as we show below, positive g)- , we have 
that C2(to, n) > below the critical dimension. 



The Fixed Point and the Upper Critical 
Dimension 



Let us have a closer look at the couplings for two- 

15(1,2) 

2 9r 



Xr- From the 



particle interaction, jg?'"^^ 
flow equation for the rescaled potential, Eq. (13 1, we ob- 
tain 

d^Xr = {d-2)Xr + 2VdXl. 

At the upper critical dimension dc = 2 there is a trans- 
critical bifurcation, such that, when the dimension d < 2, 
there is an unstable flxed point at A7- = (recall that r 
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flows in the negative direction) and a stable one at 



(17) 



For d> 2 the stabiHty of the fixed points is interchanged 
and at d = 2 they merge to one, marginally stable fixed 
point. 

Similar behavior is observed for the other rescaled co- 
efficients qi™'"-*. Once the coefficients which pre- 



cede Or ' m tne order oi tne calculation indicated m 
Fig. [2] have relaxed to their fixed point ), their 

flow is described by an equation of the form 



in the order of the calculation indicated in 



m.n) 



Ci(to,7i) +C2(m,n) • , 



with strictly positive C2{n,m) when d < 2. Hence (jtI™'"-* 
approaches the flnite and stable flxed point ^*(™'") 
— ci(to, n)/c2(TO, n). By induction this holds for all non- 
vanishing coefficients ^*(™'"). 

Thus, below the critical dimension, the flow drives the 
rescaled potential to a flxed-point potential Ur — > u*, 
which can be represented in the form 



E 

m>l,n>2.m<n 



---*-(m,n) -m n 



X X 



(18) 



In contrast, above the critical dimension, Ur tends to 
zero. In this case, we consider the dimensionful potential 



Uk instead (see Section IV I . The critical dimension dc = 
2, where both potentials, Ur and Uk tend to zero along 
the flow, is treated separately at the end of this section. 



C. The One-Dimensional Case 

Simple scaling arguments (see e.g. [12]) already indi- 
cate that the density will behave as 



p ^ Ar 



(19) 



in the long-time limit, when the dimension rf = 1, for 
some amplitude A: The density p corresponds to the fleld 
-0, such that under renormalization it scales as p = kp, 
with the "dimensionless" density Pj_see Eq. ( [TT| , whereas 
time scales a,s t = k~'^t, see Eq. (10 1. In the following. 



the most difficult task is to estimate the amplitude A. 

We deflne the rescaled non-equilibrium force by 
i^fc(V') := 9^J7fe(?/', V')|^=o its dimensionless counter- 
part by frix) ■= t^x'"i"(X; x)lx=o- Just as the rescaled 
potential Ut flows to u* , the renormalization group flow 
drives to its fixed-point value /*, which according to 



Eq. ( 18 ) may be written as 



n>2 



The kinetic equation becomes 
dtp=-\imk^fr{k-^p)^' 



*(l,n). 



limfcV*(fc"V), (20) 
k—^0 



where the second equality is valid to lowest order in k. 
The limit must not depend on k, since, once the recipro- 
cal scale k~^ is much larger than the correlation length, 
the right hand side of the equation should have converged 
well. Hence, at the fixed point we will have /*(x) ~ cx'^, 
when X is large, for some universal factor c. This implies 



that the non-equilibrium force F{p) 
kinetic equation (15 1 becomes 



and that the 



dtp = ~cp^ 



such that we indeed recover the decay law, Eq. ( 19 1 , with 
A^{2c)-i. 

Determining the factor c is tantamount to calculating 
f*{x) for large values of x- This in turn affords a good 
knowledge of the fixed-point potential u*{xjX)- Typi- 
cally, the goal of the numerical calculations is to extract 
critical exponents by considering the fiow in the region 
around the fixed point. In this case, to obtain a satisfac- 
tory result, it is often sufficient to perform a series expan- 
sion of the Wetterich equation to the first few orders in 
X and X and then to consider the flow of the coefficients 
~<jn,n) problem this clearly will not suffice, since 

the lower order coefficients only describe the behavior of 
the force /* around the origin but not for large x- 

We have exploited the special simplifications in the 
flow for the coagulations process to calculate a large num- 
ber of flxed-point coefficients g*!™'"). The equations were 
solved exactly (yet of course within the truncation of 
Eq. (|8|) employing computer algebra software. We were 
thus able to extract the flrst 125 coefficients 
power series of f*. The behavior of f*{x) for large x was 
evaluated in a double logarithmic plot, cf. Fig. [3j Since 
the power series has a flnite radius of convergence we en- 
hanced the result by employing Fade extrapolation |43j . 
For large values of x, the terms in the expansion indeed 
add up a to a power law of the order x'^- We flnd that 
approximately 



/(X)^4.2x=^ 



p(t) - 0.35 r 



1 

27re 



As compared to the perturbative result A = 
^'"^1^^"^ « 0.22 of [B] (with e = dc-d = 1), this is much 
closer to the exact decay amplitude ^== « 0.40 |44fl49] . 



D. Generalization to d < dc 

Formally one can extend the above approach to "di- 
mensions" d below the critical dimension dc = 2 and 
calculate the corresponding amplitude Ad in the long- 
time scaling of the density p ~ Ad ■ t~ 2 (see Figure Wk . 
In complete analogy to the previous section, we flnd that 
for large values of the fleld x, the fixed-point result of the 
non-equilibrium force must scale as /*(x) ~ c^X^^ with 
the dimension-dependent but otherwise universal factor 
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dimension d 



Figure 3. Double logarithmic plot of the rescaled non- 
equilibrium force /t = /* (solid red line) at the fixed point 
for a one-dimensional system. It was obtained by calculating 
its expansion up to order 125 in x within the local potential 
approximation. For small x, f*{x) = f X^ + irX"^ + - ■ • is dom- 
inated by the power law ^x^ (As-t dashed line). When x is 
increased, one reaches a regime where /* is described well by 
4.2 x"^ (steep dashed line), before the Fade approximation — 
utilized to extend the regime of convergence — breaks down. 
As expected, cf. Eq. ( |20[ ), for large enough x, in addition to 
the third order term ^ ■ x^y rest of the expansion of /* 
combines to another term scaling as x^- 



Figure 4. Results for the amplitude Ad of the long-time decay 

p ~ AdP~^ (d < 2). The circles are our estimates from the 
non-perturbative renormalization group calculations within 
the local potential approximation. They are compared to the 
exact result = —i^ in one dimension [44 and to the per- 

1 , 21n(87r)-5 



turbative result ^ — h 



1 

^'•"'1;'-^ of [in] (solid line). The latter 
IS an expansion m the deviation e from the upper critical 
dimension, e = dc — d. It is doubtful whether e = 1 (corre- 
sponding to one dimension) can be considered small. Indeed, 
despite the relatively crude truncation, the non-perturbative 
approach provides a substantially better result. 



Cd- Hence for the kinetic equation we have 



E. Treatment for the Critical Dimension 



dtp = -F{p) 



-CdP 



As for the one-dimensional case, the factor Cd is gleaned 
from a sufficient number of coefficients in the expansion 
of the rescaled non-equilibrium force friXiX) — f*{XiX) 
at the fixed point. One observes that the Pade approxi- 
mation works the better, i.e. converges for larger values 
of Xi the closer one approaches the critical dimension. 
Indeed, from perturbation theory one expects that near 
the critical dimension only and g^.^ impor- 

tant, so that the approximation should converge quickly. 
Performing the limit dc^d = e ^ 0, we can make contact 
with a result from perturbation theory |16| . To this end, 
we assume that to lowest order 



f{x)-CdX'- 



i.e. the constant Cd is, to good approximation, equal to 
the coupling A* — ig**^^'^^. This is plausible since the 
exponent l^f ~ 



sa 2. Our assumption indeed allows to 
derive from Eq. (17 1 the relation Cd = 27re + O(e^), and 
thus we recover the result from perturbation theory Ad = 
7^ to leading order in e. 



At the critical dimension dc = 2 the couplings 
5t^'^\ 5r^'^'' behave as ^ when r —oo. Since the Feyn- 

man diagrams which determine the flow of g^™'"^ involve 
n of these "elementary" couplings, the other coefficients 
go to zero as Therefore, the potential vanishes along 
the renormalization group flow, Ur — >■ 0, we cannot take 
the limit of Eq. ( 20 1 and a straightforward application 



of the analysis of the previous sections to determine the 
long-time behavior of the density is not possible. In- 
stead, we start, at finite renormalization time r, with 
the dimensionless equation 



dix = -[x + fr{x)] {i^kh). 



(21) 



where the constant term x stems from the dimension- 



less cutoff function (12 1 at vanishing rescaled momen- 
tum q = {limg^o (pr{q^) — 1). In this equation the 
long-range fluctuations have not yet been integrated out 
completely. At finite t the results roughly correspond 
to those of a system of finite size, with edge length k^^ 
(recall that k = Aexp(r)). 

In the previous subsection we saw that, as the critical 
dimension is approached, d — > 2, the decay amplitude Ad 
is determined by the lowest order coefficient A^. = ^5t^'^'* 
(up to corrections in the difference e = 2 — d). Therefore, 
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let us assume that we may set 

T 

ignoring higher order coefficients (for a more rigorous 
analysis by means of the perturbative renormalization 
group, we refer to |16|). A constant initial density po im- 
plies a diverging dimensionless initial density ^ . Solving 



Eq. (21) for this initial condition, we obtain 



X{i) = 



Pqt 



as long as X is large enough so that the term — % in 
Eq. (21 1 can be disregarded. (It destroys the algebraic 
behavior for large times, where x{t) ^ exp(— t) decays 
exponentially, as one would expect for a finite-size sys- 
tem.) Thus, after an initial transient time of the order 
rfc^, which goes to zero exponentially fast as r — — oo 
and can therefore be neglected, we have y(<) = — tt-t. 

Inserting x{i) = ^ 
ip = p) and ' 

suit [isi 



and t 
\n(k/K) ' 



k^t (see Eqs. (10 



Mi) 

2 ' 



111 with 



we recover the re- 



Ini 
4^ ' 



for large times t. 




time 




Figure 5. 
grams for 



In (a) and (b) the one-loop Feynman dia- 
the calculation of the flow of the decay rate 



2yfc 2*- fc (O,0;O,0;O,0) 

Diagram (a) determines the flow of the (l,2)-vertex func 



,(27r)'''(y"r)-' are shown. 



tion r 

r(2.2) 



(1,2) 
fe (O,0;p,i 



and involves a (2,2)-vertex function 

i , . , , , , whose flow (deduced from dia- 

gram (b)) provides a closed and exact formula for the (2,2)- 
vertex function. Diagram (c) determines the flow of the co- 
eflicient g\' . The three internal propagators imply a factor 



(j^) for the behavior oi g\, in the limit of small k. Inte- 
gration over momentum and frequency space in the Wetterich 
equation abates this singularity by a factor W^^"^ . Hence, when 
2 < d < 4 the coefHcient g^^'^'' diverges as k'^~^ . Above four 
dimensions it approaches a finite value as fc — >■ 0. 



IV. BEHAVIOR ABOVE THE CRITICAL 
DIMENSION 



We have seen that below the critical dimension two 
the renormalization group flow drives the (dimensionful) 
renormalized decay rate Afc := \9k '^^ to zero. In or- 
der to treat this singularity, we introduced the dimen- 
sionless couplings — fc-2+d(n-i)^(™,n)^ which tend 
to a finite, universal value the infrared cut- 
off scale k becomes small. Without tuning of parame- 
ters, the couplings flow to a fixed point, which implies 
anomalous long-time behavior p ~ Adt~ ' with a univer- 
sal amplitude Ad- We shall show in the following that in 
contrast to this, above the critical dimension the renor- 
malized decay rate Aj, attains a finite value as the scale 
k goes to zero. Thus, above the critical dimension the 



long-time behavior obeys the LMA dtp 



with 



non-universal macroscopic decay rate p. :— \k=o and at 
long times we recover the "classical" scaling p ^ p^^t^^ . 



A. Derivation of the Macroscopic Decay Rate 

The renormalized reaction rate can be obtained from 
the identity 



1, 



where y = = {2-kY5{0) and T = Jdt ^ 2t:6{0) de- 
note the asymptotically large volumes in space and time 
which are summed and integrated over, respectively (in 
the final result they drop out). 

Let us first restrict to the case where there is only lo- 
cal interaction between the particles, A(y — x) = (5x,yA, 
so that they can be regarded as extending over one site 



of the lattice. To calculate F 



(1.2) 
k (O,0;0, 0:0,0) 



we consider the 

.(1,2) 



flow of the more general vertex function F, ,„ „ ^ ■ 

^ k ( U,0; — p,— a;;p,tJ ) 

It is determined by the one-loop diagram pta) and apart 
from a (l,2)-vertex contains also a (2,2)-vertex. The flow 
of the (2,2)-vertex function T , J , , , , . 

V ' ^ (p ,0; , — p , — cj ,p,cl;; — p,— cj) 

follows from diagram [Sjjfo). It is self-contained in the 
sense that it only depends on the (2,2)-vertex function 
itself. Since the internal momenta and frequencies of the 
one-loop diagram do not depend on the external ones. 



(2 2) 

we have that the identity FJ, 



-(2,2) 
k (O,0;0,0;0, 0:0,0)' 

to all scales k. 

T-(l,2) 



which holds at k 

By the same 

_ -r(2,2) 



— p', — C(;':p,a;: — p,— oj) 

= cx), is bequeathed 
token, the identity 
(2ttI — which 



A: (O,0;p,ti);-p,-w) fc (0,0:0, 0;p,w;-p,-tj) 2 

evidently holds for the microscopic action ([!]), remains 
valid along the renormalization group flow. Thus, evalu- 
ating diagram [5]J6), we obtain the flow for the renormal- 
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ized decay rate 

dk^k — — 29*; 



A? 



(i?fc(q) + f(q)r 



2X1 



4e(fc2-6(q)) 



/c3 



Therefore, in the case of one-site objects the exact solu- 
tion to the macroscopic decay rate reads 



111 
^ Aq A 



1 

q e(q) 



This equation connects the microscopic reaction rate A 
with its macroscopic counterpart /z by a term which de- 
pends on the structure of the lattice. The solution is fi- 
nite on condition that the dimension d > 2. In contrast, 
when d < 2 the integral diverges to infinity, indicating 
an anomalously slow decay, which was discussed in the 
previous section. For a cubic lattice with lattice spacing 
a = 1 and diffusion constant D = 1, by numerical in- 
tegration we find that = X'^ + 0.252731009858(3). 
This value is corroborated by our numerical simulations, 
where we have considered the long-time decay of the den- 
sity p ~ ^-^-i [T7]. 

We now proceed to derive the flow equation to the 
renormalized reaction rate A^ for general reaction ker- 
nels A(z), whose interaction may extend over several 
sites. Similar as for one-site objects, we have the identity 



r 



(1,2) 



r 



(2,2) 



(271-)''+l rp, 

again it suffices to calculate the (2,2)-vertex function 
ri ; , , , s, with the important difference, 

however, that in general this vertex function depends ex- 
plicitly on both of the momenta p' and p. 



-(2,2) 

k {p' ,Lo' -. — p' .,—uj' -jp.uj-j—p. — uj) 



AVT 

(271)4^+4 



A(p',p), 



Notice that there is no dependence on the frequencies w' 
and Lo. We obtain 

e(fc2 - e(q)) 



9feAfe(p',p) = 2 / Afc(p',q)Afc(q,p; 



fc3 



Equivalently, in position space 



afcAfc(x,y) = 2/ Afe(x,q)Afe(q,y) 

Jq 

which may be rewritten as 



9(P - e(q)) 

fc3 



, (22) 



9fcAfe(x,y) = ("Pi o Afc)(z,y) 



p X! o Afc)(x,z)Afc(z,y) , 



where the projections Vi and 1^2 are defined by 
(7'ioAfc)(z,y) = / exp(zq-z)Afc(q,y)e(fc2-e(q)) 



(T'zo Afe)(x,z) = I exp(zq-z)Afe(x,q)e(fc2-e(q)). 

iq 

For numerical calculations it is an important simplifica- 
tion that the support of Afe(x, y) is contained in S x S, 
where S is the support of A(x). 

As an example for extended objects, consider balls of 
radius ^ in three-dimensional continuum space. The 
reaction kernel A(z) = A8(i? — z), where z is the dis- 
tance between the centers of the particles. Despite the 
conservation of the support, Eq. (22) is still difficult to 
solve for general A. However, a particularly interest- 
ing case are instantaneous reactions, obtained by tak- 
ing the limit A -> cxo. This corresponds to the classical 
problem treated by Smoluchowski |10| , who argued that 
/i = AttRD (where D denotes the diffusion constant), 
by a heuristic approach, which was later rigorously con- 
firmed by Doi [8J. Notice that the problem is equivalent 
to A(z) = XS{R — z), where the kernel is nonzero only 
on the surface of a sphere, in the limit of infinitely fast 
reaction, A — >■ c». We now derive a simplified exact fiow 
equation for this and similar cases. 

Suppose that the kernel Afe(x) 
X)y Afc(y, x) can be described by only one degree of free- 
dom, i.e. on its support S, the reaction kernel Afc(x) is 
a constant (for fixed scale k) and all the elements of the 
support are equivalent, in the sense that every element of 
the support can be mapped onto every other element un- 
der rotations and refiections that conserve the action S. 

According to Eq. (p3| 



Ey A/c(x,y) 



afcAfe(x) = ■^'^Xki:x.,z)'^{'Pi o Afe)(z,y) 



fc3 



5]Afc(x,z)(7'oAfe)(z). 



(23) 



with the projection 

(V o Afc)(z) = f cxp(iq • z)Afc(q)e(fc2 - e(q)) . (24) 

J q 



Since (V o Xk){z) is a constant on S, the sum in Eq. (23 1 
is trivial and we obtain 

4A.W. ^^-M';:,°^''W , (25) 



Alternatively this can be written as 



dkXkip) = 2 Afe(p-q)Afe(q; 



e(fc2-e(q)) 



(26) 



As elaborated in the Appendix, from Eq. ( 25 1 one ob- 



tains an analytic solution for the kernel A(z) = XS{R— z) 
in continuum space. In particular, we recover Smolu- 
chowski's classical result /i = AttDR for the macroscopic 
decay rate of spheres of radius R, diffusing with diffu- 
sion constant D and coagulating upon contact. In the 
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stochastic simulations we found it preferable to work with 
objects defined on a lattice rather than in continuous 
space. In Fig. |7]we give two examples of extended object 
defined on a lattice, whose exact macroscopic decay rate 
can be derived from Eq. ( 26 ) . 



B. Universal Correction to the Law of Mass Action 

To study corrections to the LMA, let us represent the 
non-equilibrium force as the limit of a power series 



1 



n>2 

exploiting the fact that F^ is analytic if k > [26] . The 

(1 2) 

lowest order coefficient = 2Afe, treated above, con- 

verges to a constant value, the macroscopic decay rate 

(12) 

5fc=o ~ Thus, assuming that the higher order terms 
^\9k^'^^ (n > 2) can be neglected, we would recover the 
LMA term F{p) = /xp^, quadratic in the density p. It is 
crucial to notice that the next to leading terms are not 
-,(1:3) ^3_L ji^„(i,4) ^4 Rather, in d dimensions 

with n > turn out to diverge 

as k goes to zero. In the following analysis, we show that 
the infinite sum of these diverging terms converges and 



simply i g^t^o P^ + h. 9k=a + • 
all coefficients g'\!''"'\ 



gives the finite contribution X]ri>^±? 'h.9k '^^^ P" 



CdP' 



(up to possible logarithmic factors in p), for some con- 
stant Crf. In three dimensions this provides a relatively 
large leading correction. 

The flow of the physically relevant couplings 5^"'^'"'' is 
calculated from diagrams with n incoming and one outgo- 
ing leg. To lowest order, their divergence stems from the 
contribution to the flow of diagrams which only contain 
(1,2)- and (2,2)-vertices and follows from power counting. 

(1 3) 

Let us exemplify this for the (l,3)-coupling g^,' , with 
the associated one- loop diagram of Fig. [sjjc). For any 
finite scale fc, also the couplings must be finite. There- 

(13). 

fore, the divergence in k of gj, ' builds up in the limit of 
small fc, where only long- wave length and short-frequency 
fluctuations, q < fc, and uj < fc^, contribute to the flow. 
Above the critical dimension, the (1,2)- and (2,2)-vertex 
functions are not divergent for fc = q = w = 0, but attain 
a finite value, which is equal to the macroscopic decay 
rate p. Hence we can take their value at k ~ q = ui ^ 0, 
if we are only interested in the strongest divergence. The 
evaluation of the diagram then yields 



dkgi'^^^ = -4l6 / Gk{<i,ujfGk{-q,-uj){~pf 



with the propagator 



(27) 



(fc2 - e(q)) e (fc2 - e(q)) + e(q) + ' 



and the macroscopic reaction decay rate p — /i(q = 0). 
Due to the derivative dk = dkRk ■ dR^i the integration is 



restricted to the domain e(q) < fc^, where the propagator 
is independent of q, i.e. 

dk =dk 

such that within the domain of integration the propaga- 
tor simplifies to 



1 



k'^ + iui 



The fact that the integral, Eq. ( p7| ), does not depend 
on the full reaction kernel p{q) already indicates that 
these divergences cannot depend on the shape and size 
of the objects: Originating in long wavelength fluctua- 
tions around q = 0, they do not resolve the details of the 
reaction kernel. 

To lowest order in fc, the dispersion relation e(q) can 
be approximated by the "continuum limit" e(q) = q^. 
Therefore, the divergences are not only unaffected by 
the shape and size of the particles, but also independent 
of the structure of the lattice; it is as if the divergent 
terms only "see" structureless point particles (for which 
/i(q) = p, independent of the momentum q) that are 
embedded in continuous space (where the dispersion is 
simply e(q) = g2). To lowest order in k the integration 
over the momentum q then yields the volume of the d- 
dimensional ball with radius k. Thus, up to some positive 
constant factor, Eq. (27 1 becomes 



3ud 



p-'k 



{l + ivY{\-w) 



-dk 



The integral, where we substituted ^ — v, contributes a 

positive factor. Thus, the coupling g^^''^'' scales as k'^~'^ 

for dimension c? < 4, diverges logarithmically at c? = 4 
f 1 3) 

(in both of these cases g\. ' diverges to positive inflnity), 
and converges to a flnite value for c? > 4. 

In summary, and generalizing to arbitrary couplings, 
the strongest divergence of a diagram is obtained by the 
following prescription. Each propagator gives a factor 

Gfc(q,w) - 

The resulting divergence is attenuated by the integration 
over the momenta q and frequencies 



d'^q - fc^ / dw 



(28) 



where in the left hand formula the exponent 3 is simply 
the dimension. Finally, the lowest order contribution in k 
of the vertices must be multiplied with the result, in par- 
ticular (1,2)- and (2,2)-vertices give rise to the constant 
factor — /i. 

It is now straightforward to determine the strongest 
divergences of the couplings g^^'"''. We need to con- 
sider the one-loop diagrams with n incoming and one 
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outgoing leg that contain only (1,2)- and (2,2)-vertices. 
Clearly, in these diagrams there are exactly n of these 
vertices, connected by n propagators, which gives a fac- 
tor (— /i)"'fc~^"'. After integrating over the momenta and 

frequencies, cf. Eq. (28 1, we have that the couplings g^'"''' 

converge to a finite value go^'""* if 2n — (d -I- 2) < and 
otherwise diverge. 



(l,n) 



9 k 



if 2n<d + 2. 



d + 2, 

(-l)"+icrf.„,M"A:''+2-2" a 2n>d+2. 



(-l)"crf,„Ai"ln(A:) if 2n : 



(29) 

with some Cd^n > 0. 

Diagrams with higher order vertices (more than two 
incoming legs) can be neglected. Suppose, for instance, 
that such a one-loop diagram contains a (1, n')-vertex 
(n' > 2), which contributes a factor /c^^+^-S" ^2n' > 
d + 2). If we replace this vertex by a string of rt' — 1 
(l,2)-vertices, connected one by one by rt' — 2 propaga- 
tors, this results in the more relevant factor of the order 
^-2(ti -2) (jQQg^ii i^jjat the dimension is larger than the 
critical dimension dc = 2). 

The fact that the couplings g^^'"'' diverge in an alter- 
nating sequence is not surprising, given that one also ob- 
tains such sequences when one takes the Taylor expansion 
of the non-analytic function a;" for a non-integer a > 
around, say, x = 1. Indeed, we are on the lookout for 
such non-analytic terms. The infinite sum of diverging 
terms can be written as 



E 

2n>d+2 



(30) 



for some scaling function fd- Since for large systems the 
non-equilibrium force must become independent of the 
system size, corresponding to the reciprocal scale k~^, 
one obtains 

fd[X) ^ CdX , 

for some constant Cd- In particular, this indicates that in 
three dimensions, as opposed to higher dimensions, the 
next to leading term to the LMA term in the force F is 
not of the order three but a non-analytic term of order |, 

F{p) = fip^ + csi^ip)^ + ... 
In fact, we cannot strictly rule out that the sum in 



Eq. (30) does depend on k. The fc-dependent can cancel 



with some other correction term (or terms) to the poten- 
tial. Actually, this can rectify the problem of the loga- 
rithmic correction in even dimensions. Assuming that 



fk (x) 



Cd ■ x^^' 



ln(x) , 



for even dimensions d, then the sum in Eq. ( 30 I , in ad- 
dition to a fc-independent term Cdipip)^'^^ ln{pip), gives 
rise to a term —2cd{pip)^'^^i'D.{k). This term can cancel 
with the term that is logarithmic in k, cf. Eq. (29 I. 



Let us finally exploit our findings to calculate the cor- 
rection term exactly from the Wetterich equation. We 
choose the ansatz 



(31) 



with the diffusion term 5f and the term Sz=i as defined 
in Subsection III Al The initial condition for the effective 
average potential reads Uk=oo{'4',''P) = A(q = O)^^-!/)^ -|- 
A(q = 0)i}ilj^, with A = A(q = 0) ^ X^x •^(^)- According 
to our above discussion, Eq. ( [3T| includes all the terms 
needed to determine the correction exactly, even after 
subsituting 

Ukiip, V) ^ mV'^V'^ + mV'V'^ , e(q) ^ ■ 

For the renormalized non-equilibrium force Fkitjj) — 
7°''^'* we then obtain 

dip 



dkFkip) 



2Vdk^+'p^p^ 
(fc2 + 2ppy 



(32) 



where Vd stands for the volume of the d-dimensional 
sphere with radius {2tt)~^. Although this equation is 
not exact, it still delivers the correct non-analytic contri- 
bution of the order p^~ , because it treats all the terms 
which give rise to it exactly. Indeed, integrating Eq. (32 1 
from any fc = A>OtoA: = yields a contribution 



MP 

2 



(where F denotes the F-function and not the average ac- 
tion) which is valid for a dimension d > 2, as long as it 
is not an even natural number. We notice that the cor- 
rection alternates its sign, from a positive contribution in 
three dimensions, to a negative in five dimensions, and so 
on. As d approaches an even number, the result diverges, 
indicating logarithmic correction terms. In four and six 
dimensions, for instance, we find that a term 



87r2 



327r3 



respectively, is added to the force F. The fc-dependent 
logarithmic term in Eq. (29) indeed cancels. We remark 



that the correction terms are also amenable to perturba- 
tive treatment |T7, "50]. 

We have run simulations for a range of different models 
(one-site objects with both finite and infinitely large re- 
action rates, and two examples of extended objects that 
react immediately on contact) for the three-dimensional 
system, where 



F{p) = pp^ + 



p-2 



2V2tt 



(33) 



The simulation results clearly corroborate our theoretical 
findings (see Fig.[6|. 
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Figure 6. Rescaled data for the universal correction to the 
non-equihbrium force F. In the stochastic simulations, F 
was determined directly by introducing homogeneous particle 
input and considering stationary states. On this double loga- 
rithmic plot, we show the rescaled data for {F{4>) — ^^P^) /m' 
for a range of models. We predict this term to be of the 
universal form / (2-\/ 27r) (solid black line), independent of 
the model, cf. Eq. ( 33 1 . The data evidently corroborate our 
theoretical results. 



V. CONCLUSION 

In this article we presented results of our study of 
the coagulation process by means of a non-perturbative 
renormalization group (NPRG). Below the critical di- 
mension the renormalization group flow drives the di- 
mensionless non-equilibrium force fr (%) to a unique fixed 
point /*(x). Within a certain approximation we can cal- 
culate f*{x) thus derive the ensuing anomalously 
slow, universal long-time density decay p{t) ^ Adt^^- 
As the dimension is lowered from 2, the universal am- 
plitude Ad becomes a non-perturbative quantity, beyond 
the reach of the perturbative approach. By considering 
the power law regime of /*(x) for large Xi the NPRG 
enables us to extract the universal amplitude Ad even 
in one dimension where the estimate for Ad is in good 
agreement with exact calculations. 

Above the critical dimension, the long-time decay is 
governed by the law of mass action (LMA), such that 
dtp = —pp^ to lowest order in p, and p{t) ~ p,^^t^^ in 
the long-time limit . NPRG provides a closed formula for 
the macroscopic decay rate p. Starting from the micro- 
scopic rate (which may be infinitely large) fluctuations 
in space and time are integrated gradually, going from 
short wavelengths and frequencies to long ones. Along 
the renormalization group flow the effective decay rate 
becomes monotonously smaller until all contributions are 
integrated and one obtains the macroscopic decay rate 
p. We work out the solution to the flow equation for the 
decay rate for a number of examples which allow for a 
particularly accurate solution. 



Furthermore, we find that there are correction terms to 
the non-equilibrium force F violating the law of mass ac- 
tion, which is generalized to dtp = —F{p) = —pp^ + . . . 
We identify a non-analytic term Cd{pp)^^ in the non- 
equilibrium force and show that it originates in long- 
range and short-frequency fluctuations. The term is uni- 
versal in the sense that the amplitude CdP^^ depends on 
the particular features of the process only through the 
non-universal rate p. The factor Cd is completely inde- 
pendent of the microscopic details of the reaction process. 
For the three-dimensional case, we have run stochastic 
simulations which clearly confirm the theoretical predic- 
tions. 

The NPRG is a versatile and powerful tool for the 
study of non-equilibrium systems and the coagulation 
process is not only a simple model for the analysis 
of more complex systems but also serves as a starting 
point for the study of more complicated theoretical mod- 
els |23| |25| l5T| [52] . We therefore expect that our results 
are relevant for a range of experimental systems and the- 
oretical models and believe they will encourage further 
work on the fundamental implications of fiuctuations on 
non-equilibrium processes. The strong impact of fluctua- 
tion on the coagulation process below the critical dimen- 
sion has already been probed by experimental studies 
on effectively one-dimensional exciton dynamics [551 [M] . 
Similar experiments on excitons which disperse in all 
three spacial directions [55 may be suitable for analy- 
sis of our predictions of a violation of the LMA in three 
dimensions. 
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Appendix: The Macroscopic Decay Rate for 
Selected Reaction Kernels 



In the following we study reaction kernels which allow 
for a particularly precise numerical solution, repeating 
for completeness the calculations given in the Supple- 
mentary Material of [T7]. Consider the reaction kernels 
whose two-dimensional versions are depicted in Fig. [7] 
We study their three-dimensional versions on a cubic 
lattice. The lattice sites that these objects consist of 
(i.e. the sites that are part of the reaction kernel) are all 
equivalent. This symmetry must hold along the renor- 
malization group flow. Furthermore, the flow equation 
conserves the support of the reaction kernel in position 
space. Thus the objects also retain their shape. Explic- 
itly, in three dimensions the renormalized reaction kernel 
of extended object 1 can be expressed as 



Afc(p) 



3 

: 2Afc ^ cos(p^) , 
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Figure 7. Two-dimensional versions of the reaction ker- 
nels of extended objects (solid red). In three dimensions, 
for extended object 1, the kernel A(z) = Xoo if z G S = 
{(±1,0,0), (0, ±1,0), (0,0, ±1)}. Otherwise it is zero. Notice 
that for instantaneous reactions, the striped square can be re- 
garded as part of extended object 1, which is then a discretiza- 
tion of the sphere. The support S (with A(z) = Aoo if z G S) 
of the three-dimensional reaction kernel of extended object 2 
is crated by the union of the set {(1, 1, 2), (1, 2, 1), (2, 1, 1)} 
with its mirror images in each octant. Since the flow con- 
serves the support of the reaction kernel, their shape re- 
mains the same. For instantaneous coagulation reactions 
we find fi'^ = 0.086064343192(3) (extended object 1) and 
/i"^ = 0.036287603611(2) (extended object 2). 



The microscopic decay rate is A :— Aoo(O) = 6Aoo (Aoo 
is the reaction rate for each site of the kernel). Solving 
Eq. (26 1 at p = 0, one obtains the macroscopic decay 
rate 




1 _ 1 _ 1 
M ~ Ao(0) ~ Aco(O) 



e(fc2-e(q)) 



/fc3 



1 

A 



E,.=iCOs(p,.) 



q36E.=iSin'(9./2) 



0.086064343192(3) 



In the last step we inserted the dispersion relation e(p) = 
4 ^j^^j^ sin^ (^) for the cubic lattice with unit lattice 
spacing. The calculation for extended object 2 is analo- 



gous. If each of the 24 sites effects a reaction with rate 
Aoo then 



- = I ±0.036287603611(2). 
/i A 



with the microscopic decay rate A — 24Aoo- The simu- 
lations confirm these results for the macroscopic decay 
rates and support our prediction on the universal correc- 
tion to the non-equilibrium force, cf. Fig. [6] 

Finally we choose the reaction kernel to be the sur- 
face of a sphere, Afc(z) = Xk5{R — z), for objects dif- 
fusing in continuous space. In the limit of an infinitely 
large microscopic reaction rate this is evidently identical 
with spheres that coagulate instantaneously on contact. 
Without loss of generality, we set the radius i? = 1 in 
the following. Using spherical coordinates, the projec- 
tion {V o Xf.) (z), see Eq. (24 1, becomes 



1 



(27r)3 



dq dl9 d<P q^ sin(i?)e^«" ''°'^'^'> e {k^ - q^) 

dr dd d^ sin(?9)e-'"f ™«('') Xk5{l-r) = 
2 f'' 1 

— / dq - sin(gz) sin(g) . 
Jo z 



= --fk(z) 



Thus from Eq. (25 I we have dk^k = 

/fe(l) 



2A2^ and 



1 

Ao 



1 



dk 



I 



1. 



1 

4n ' 



Aoo -^0 ^'^ Ao 
The macroscopic decay rate becomes 
1 _ 

where A — 47rAoo- Therefore, for instantaneous reactions 
H = AttDR , 

with diffusion constant D and radius R. This confirms 
Smoluchowski's result [TU], proved to be exact by Doi [5]. 
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